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. Abstract 

^ We investigate tiie transition plienomenon of the universe between a phantom and a non-phantom 

I phases. Particular attention is devoted to the case in which the cosmological constant depends on 

I time and is proportional to the square of the Hubble parameter. Inhomogeneous equations of state 

' are used and the equation of motion is solved. We find that, depending on the choice of the input 



o 

u 



parameters, the universe can transit from the non-phantom to the phantom phase leading to the 
appearance of singularities. In particular, we find that the phantom universe ends in the singularity 
of type III, unlike the case without variable cosmological constant in which the phantom phase 
ends exclusively in the big rip (singularity of type I). The Cardy-Verlinde formula is also introduced 
O . for inhomogeneous equation of state and we find that its equivalence with the total entropy of the 

. universe, coming from the Friedmann equations, occurs only for special choice of the input parameter 

■ , . 

m at the present time. 
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. 1 Introduction 
rn 

I Accelerated expansion seems to play an important role in the dynamical history of the universe. There 

is a firm belief, at present time, that universe passed through inflationary phase at early times and there 
are growing evidences that it is accelerating at present. The study of large scale structures indicates 
' that the universe is almost spatially flat and that dark energy accounts for about 70 percent of the total 

^ I energy content [T1[11I31I1]- Moreover, dark energy is believed to be the responsible of the acceleration of 

our expanding universe. This sort of fluid violates the strong version of the energy conditions. Besides, 
when the null version is also violated, the fluid is called phantom and then the universe may present 
future singularities at finite time [3 [?]• In the phantom phase, the energy density grows whereas it 
decreases in a non-phantom one. However, we know little about the nature of dark energy in general and 
of phantom fluid in particular, except for their negative pressure. Therefore, a large effort has been spent 
in recent years to explain this mystery [8j [9l 1101 llli 1121 113) . Specially, this transition phenomenon known 
as quintom scenario is proposed by B. Feng and collaborators in [T4j, and they found that it gives rise to 
the equation of state larger than —1 in the past and less than —1 today, satisfying current observations. 
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For understanding the possible connections among the dark energy models, it is useful to study the cosmic 
duality. Then, it was studied the duality in two-field quintom models of dark energy and it has been found 
that an expanding universe dominated by quintom-A field is dual to a contracting universe with quintom- 
B field jl5) . Recently, Yi-Fu. Cai and collaborators wrote a paper which introduced the experimental 
developments on finding the transition between quintessence and phantom phases and various theoretical 
realizations of such a scenario, see [TS]. On the other hand, the knowledge of some properties of the 
universe in the phantom phase, such as the singularities, motivated various investigations with the aim 
of dealing with them. A. B. Batista and collaborators 117^ investigated the effects of particle production 
when a massless minimally coupled scalar field is present in spacetimes where a; is a constant. To do 
so they used a state for which Bunch and Davies [18j had previously computed the stress-energy tensor. 
They found that the energy density of the created particles never dominates over the phantom energy 
density. In the same way, quantum effects near the big rip is studied in I19j where they used the n-wave 
regularization for calculating the energy density of particle creation and found that, in this case, it tends 
to infinity when the big rip is approached and becomes the dominant component of the universe. This 
means that the big rip can be avoided by a scalar massless field. Pavlov [20 computed both the number 
density of created particles and the stress-energy tensor for a conformally coupled massive scalar field 
for the case in which u — —5/3. It was found that quantum effects are not important for masses much 
smaller than the Planck mass and times which are early enough that the time until the Big Rip occurs 
is greater than the Planck time. J. D. Bates and P. R. Anderson [21] used a background field approach 
in which the energy densities of the quantized fields are computed in the background spacetime which 
contains the Big Rip singularity. They found that for fields in realistic states for which the energy density 
of the quantized fields is small compared to that of the phantom energy density at early times, and for 
spacetimes with realistic values of a;, there is no evidence that quantum effects become large enough to 
significantly affect the expansion of the spacetime until the spacetime curvature is of the order of the 
Planck scale or larger, at which point the semi classical approximation breaks down. Also in order to 
deal with singularity problem, it has been considered by Yi-Fu Cai and collaborators [22 the cosmology 
of the Higgs sector of the Lee- Wick Standard Model, an alternative to supersymmetry to solving the 
hierarchy problem. They found that homogeneous and isotropic solutions are non-singular and then, the 
Lee- Wick model can provide a possible solution of the cosmological singularity problem. 

One of the phenomenological ways to explain the dark energy problem is assuming a variable cos- 
mological constant. The cosmological constant A is pretty compatible with observation and effort is 
currently devoted to the investigation of the theoretical foundations of a variable cosmological constant 
and its model properties [23l [241 1251 [Ml- In this respect, various ansatz have been used. For example 
in [27], S. Carneiro and collaborators considered a cosmological constant A cx iJ, and studied a possible 
way to distinguish the validity of this scenario from the standard one. Note that a model A cx has 
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earlier been proposed in [28], requiring that the cosmic density p equal to the Einstein-de Sitter critical 
density pc, leading to a close universe without singularity, horizon, entropy and monopole problems. 
Furthermore, a large number of phenomenological A models have been constructed in order to describe 
the dynamics of the universe |29j . The case which attracts our attention in this paper is A cx and has 
been studied in several other works for other purposes, but always with the aim of clarifying some grey 
areas in cosmology [301 ISIl [321 123] ■ In this paper we propose to use A cx , for analysing the evolution 
of the universe, being in the phantom or non-phantom phase. The phase transition of the universe will be 
studied considering a model in which dark energy is described by some rather complicated ideal fluid with 
an unusual equation of state (EoS) which will be chosen to be an inhomogeneous one. It is important to 
emphasize that this inhomogeneous EoS corresponds to a pure dark energy models. This kind of models 
may reproduce late-time acceleration, but it is not easy to construct a model that keeps untouched the 
radiation and matter dominated epochs. However, it is easy to introduce into our considerations ordinary 
matter and radiation but in that case they only appear suddenly at some point. 

Note that this kind of study has been also considered in [34] where the cosmological constant is a 
linear function of time. Here, we use the ansatz A cx with the inhomogeneous EoS. We find that, 
depending one the choice of the input parameters of the EoS considered, the universe may transit from 
a non-phantom to a phantom phase leading to finite time singularities. Another interesting point of our 
result is that, in contrast to the model without variable cosmological constant in which the phantom 
universe ends with the singularity of type I (big rip), the phantom universe in this case ends with the 
singularity of type HI. 

Another point we address here is the Cardy-Verlinde (CV) formula coming from inhomogeneous 
EoS. Verlinde 35 made an interesting proposal that Cardy formula |36j in two-dimensional conformal 
field theory can be generalized to arbitrary spacetime dimensions. Verlinde further proposed that a 
closed universe has subextensive (Casimir) contribution to its energy and entropy with the Casimir 
energy conjectured to be bounded from above by the Bekcnstein-Hawking energy and as consequence, 
one obtains a very deep relation between gravity and thermodynamics 371 . Within the context of the 
radiation dominated universe, such bound on the Casimir energy is shown to lead to the Hubble and 
the Bekenstein entropy bounds respectively for the strongly and the weakly self-gravitating universes. 
The generalized entropy formula, called the CV formula, is further shown to coincide with the total 
entropy of the universe coming from the Fricdmann equations. These results were later generalized 
[551 15^ HOI WU [m [121 [m US] . Our goal here is to analyse the equivalence between the C V formula and 
the total entropy coming from Friedmann equations assuming that the universe is conformally invariant. 
With this, we find that for the inhomogeneous EoS, the generalized entropy of the universe reduces to the 
CV formula with a special choice of the input parameter m and this does not correspond to a radiative 
universe as in the case of homogeneous EoS. Note also that this equivalence occurs exclusively at the 



3 



presents time . 

The paper presents two sections, the first showing the inhoniogeneous EoS with which the sohition 
for the Hubble parameter is found and some discussions on the input constants are put forward, allowing 
to a whole analysis of the transition phenomenon. The second section shows a brief concept on the 
CV formula with homogeneous EoS and latter, a complete analysis of the equivalence between the CV 
formula and the Friedmann equations with inhomogeneous EoS. Finally, we present our conclusions and 
perspectives. 



2 Solving the inhomogeneous equation of state 

Les us consider the universe driven by an ideal fluid (dark energy) with the inhomogeneous equation 
of state [46] 

p = L0{t)p + A{t) , (1) 

where uj{t) and A{t) depend on the time and p and p are respectively the energy density and the pressure 
of the fluid. This equation, for the case A{t) = and uj{t) as affine function of time, has been studied in 
[46l|47]. Moreover, the case A{t) 7^ as a affine function of time has been examined in [34] . 
The equation of energy conservation and the Friedmann equations are respectively 

p + 3H{p + p) = , (2) 
^H' = P , (3) 
-1 (2H + 3H^) = -p , (4) 

where — SttG, with G the gravitational constant and H = ^, the Hubble parameter where a{t) is the 
scale factor. 

Using ([!]) and ([3]), equation ([2]) is rewritten as 

p+V^K[l+uj{t)]p^/^ + V3Kp^/^A{t) = . (5) 

From now on, we suppose that uj(t) depends linearly on time and the cosmological constant A(t) is 
proportional to the square of the Hubble parameter, that is 

uj{t) = at + l3 , (6) 

= ^P{t) . (7) 
Taking into account ^ and (O, equation (O becomes 

p+ (At + B)p3/2 ^ , A = aKV3 , B = kV3(i + /3 . (8) 
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The solution ol this equation is 

p{t) = 2 , 9 

(Ai2 + 2Bt - 2Cf 

and the Hubble parameter and its rate behave as 

Hit) = ^ , (10) 

V3 (At^ + 2Bt - 2C) 



Hit) - ^.Z'^Z ^ (11) 



8k iAt + B) 
iAt^ + 2Bt - 2C)'^ 

where C is an integration constant. 

As we are dealing with an expanding universe, we need an increasing scale factor, that is a > 0. We 

know that a = Ha, then for an expanding universe the Hubble parameter has also to be positive. Solving 

a = Ha, one obtains 



a{t) = exp Hit)dt^ 



( Ak q(t) \ , , { At + B \ 

The positivity of a depends on the sign of H{t). Note that the expression At^ + 2Bt — 2C vanishes for 



ii.2 = — (S ± vB^ + 2AC) /A. The universe expands when the Hubble parameter is positive and one 
has a phantom fluid when the weak version of the energy conditions is violated, that is when p + p < 0. 
Combining equations ^ and (jH), one obtains 

P + p = -^H , (13) 

and it turns out that the phantom phase (respectively the non-phantom one) is obtained when H > 
(H < 0). Two situations are important for a whole analysis: when the constant A is positive or negative. 
• Analysis for the case A > (a > 0) 

In this case, a simple study of the Hubble parameter sign shows that one has an expanding universe 
for ti < t < t2 and a contracting one when t <ti and t > 12- On the other hand, it is easy to see that H 
vanishes for <3 = —B/A. Then, the universe is in the phantom phase (ij > 0) when t > and the energy 
density grows; for the non-phantom phase (i? < 0) t < t^, the energy density decreases. Consequently, 
the accelerated expanding universe begins with a non-phantom phase and enters in the phantom one at 
the transition time ttr = ts, the time at which the Hubble parameter and the energy are 

AkA 16v42 ^ ^ 

Htr = —F= , Ptr = 7 • (14) 

V3iB^ + 2AC) (B2 + 2AC)' 

In the non-phantom case, the energy density decreases and tends to ptr as the time goes to ttr- The simul- 
taneous divergence of pit) and Hit) appears at ti and t2. However, a(t2) = exp 
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which is finite. In fact, the phantom universe ends with a future singularity, in this case, the singularity 
is of type III (for a classification of future singularities see |5]) since the energy density and the pressure 
at this time are divergent. The graph of H{t) is shown in the left panel of Fig. 1. 
• Analysis for the case A < (a < 0) 

For t < ti and t > t2, the universe expands whereas it contracts for ti < t < t2. The first derivative of 
the Hubble parameter is positive (respectively negative) for t < (t > t^) . The expanding universe begins 
with a phantom phase which ends at ti and enters in a non-phantom phase at t2. Here the transition 
is not instantaneous, it is the contracting phase of the universe. At ii, the energy density p(ti) and 
the pressure p{ti) diverge. However, at ti, the scale factor is finite, a(ti) — exp 

Then, the phantom phase ends with the singularity of type HI. In the non-phantom phase, the energy 
density decreases and goes to zero as t — > oo. The graph of the Hubble parameter, H(t), versus time t 
is shown in the right panel of Fig. 1. 



3 CV formula from inhomogeneous EoS fluid 

This section is devoted to the application of CV formula to ideal fluids. In a flrst step, let us tell briefly 
introduce CV formula for homogeneous EoS. We consider a (n -I- l)-dimensional space time described by 
the FRW metric, written in comoving coordinates as 

where k = —1,0, +1 for an open, flat, or closed spatial Universe respectively, and dftf^_i is the metric of 
an n — 1 sphere. Then, by inserting the metric (|15p in the Einstein equations the Friedmann equations 
are derived, 

H^^-f TTP--, H = -ip + p) + -. 16 

n{n — 1) 

The total energy E of the universe is £^ = pV, with V its total volume. Since the Casimir energy may 
be include in the total energy, we consider a closed universe, k — 1. For the homogeneous EoS, p — ujp, 
with uj a constant. Then the conservation law for energy has the form 

p + nH{l+oj)p = , (17) 

which reduces to ([2]) for n = 3 and whose solution depends on the scale factor as 

pcxa-"(i+'^) . (18) 

The total energy of the universe can be written as the sum of an extensive part Ee and a subextensive 
part Ec, called the Casimir energy, and it takes the form: 

E{S, V) = EsiS, V) + ^EciS, V) . (19) 
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Under the transformations 5* — > ^5 and V with a constant ^, the extensive and the subextensive 

parts of the total energy respectively scale as (SU [S^ 



Ee{^S,^V)^^Ee{S,V), Ec{^S,^V)^e--EciS,V) , (20) 
and therefore, we have for the total energy 

E{^S,^V)^^Ee{S,V) + le--Ec{S,V)- (21) 
Taking the derivative of ([2T|) with respect to ^ and letting ^ = 1, one obtains 

<S).-(i).,--G4)-- 

Assuming that the universe satisfies the first law of thermodynamics dE = TdS — pdV, we have the 
thermodynamics relations g — —p and {^)y ~ T. Using these thermodynamics relations and Eq. 
(IT9l) . one can put Eq. (|22|) into the following form for the Casimir energy, as the violation of the Euler 
identity 

Ec = n{E+pV -TS) . (23) 

Since the total energy behaves as £' ~ a~"'^ and by Eq. (flQl) . the Casimir energy also goes as Ec ^ 
a~"'^. The FRW Universe expands adiabatically (dS = 0) so the products Eca"'^ and Eeo^'^ should be 
independent of the volume V , and be just a function of the entropy. Then, by the rescaling properties 
(|20p . the extensive and subextensive parts of the total energy can be written as functions of the entropy 
only [37], 

where /i and are undetermined constants and Att and 2tt are used for convenience. From these expressions 
for Ee and Ec, one obtains the following expression for the entropy of the universe: 



27ra 



JIUJ 



[^^VEc{2E~Ec)] . (25) 



This result, obtained in [37], reduces to the CV formula when the universe is radiation dominated, 
Lu = 1/n, that is 

" ^™ ^Ec{2E-Ec). (26) 



Let us now look to the case of the inhomogeneous EoS that we used in the precedent section and analyse 
the relationship between the CV formula and the entropy of the universe. As has been done in |53j . we 
assume an EoS expressed as a function of the scale factor and described by 

p = w(a)p + j(a) ■ (27) 
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Introducing (|27p in the energy conservation equation ()17|) . one obtains 

n{l+u}{a)) j{a) 

P (a) H PW = , (28) 

a a 

where the prime denotes the derivative with respect to the scale factor and we took t = t{a). The general 
solution of (l28l) is 



p{a) 



Q ~ n 



, with F(a) = 



(29) 



where Q is an integration constant. In this analysis, making use of ([7]) and taking into account the 
derivative with respect to the scale factor, ([1]) can be written as 



p\a) + ^ —p{a) = 



uj(a) = — h oj(a) 



Identifying with d^H]), one gets 

j(a) = 0, F(a) = 
from which one obtains the energy density as 

p{a) = Qe-^(") . 
On the other hand, using ([T^ . one can write 



y/B^ + 2AC 
K\/3 



tan 



In (^aV3(^'+2'4C)/(4«)^ 



Making use of (pO|. (|3Tt . ([33)) and (|32)) . the energy density is written as 

p{a) = Q 



cos 



In (a) 



(30) 



(31) 



(32) 



(33) 



(34) 



Note here that only for some special conditions of the functions uj{a) and j{a) CV formula (l26l) can be 
recovered . Let us assume that the present time is = and then analyse the equivalence between CV 
formula and the total entropy of the universe at this moment. Note that as t — ^ Iq, u){a) — > rn — 1, with 

,2 

m — 1 + f3 + Then, F{a) nm\n (a) and p oc a~"™. Hence, the total energy in the volume V — 
behaves as i? = pV cx a"(^^™\ which is the same behaviour for the extensive and subextensive energy 
through ([23|l and (|19p . If we assume the conformal invariance, the products i?£;a"*^™^^^ and iJca"*-™^^-* 
do not depend on the volume and are only functions of entropy. Then, we have for the extensive and 
subextensive energy. 



from which we determine the entropy as 

"27rna"("-i) 



Er 



27ra"("-i) 



S 



/pv 



^Ec{2E - Ec) 



(35) 



(36) 



Then, for m — the CV formula is recovered. However, for any m ^ ^2±i CV formula can not be 
reproduced. Note in this case that the universe does not correspond to a radiative one. 
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4 Conclusion 



We studied the transition of the universe between a phantom and a non-phantom phases. Note that 
in the non-phantom phase, the energy density decreases while in the phantom one, it grows leading to 
singularities. We focused our attention on the variable cosmological constant which has been introduced 
in the EoS, which becomes inhomogeneous. Then, we solved the equation of motion which led to the 
explicit expression of the energy density and consequently to that of the Hubble parameter. The first 
derivative of the Hubble parameter played a crucial role in this analysis since it allowed us to know 
which time interval corresponds to the phantom or non-phantom universe. The input parameter A also 
appeared to be an important one. 

In the first part, with the EoS considered and the ansatz A oc H'^{t), two important cases have been 
found. For A > 0,we saw that the universe evolves from a non-phantom phase, ti < t < ttr, to a. phantom 
one ttr < t < t2. In the non-phantom phase, the energy density descreses and goes to the energy density 
at the transition time while in the phantom phase, the energy density and the pressure grow and go 
to infinity at the finite time t2- At the same time the scale factor remains finite and we conclude that 
the universe ends with the singularity of type HI. For A < 0, we saw that the universe begins with a 
phantom phase which ends with the singularity of type HI at ti, enters in the non- phantom at t2 where 
the energy density descreses end goes to zero as t — > oo. Here, the transition from the phantom phase 
to the non-phantom one is not instantaneous; it is a time interval corresponding to the contracting phase 
of the universe. 

However, in the case in which the phantom phase precedes the contracting phase; of the universe, it 
would be interesting to study the possible avoidance of this type of singularity introducing either the 
viscosity term in the cosmic fluid or taking into account quantum effects. On the other hand, the same 
analysis can be done with the ansatz that the variable cosmological constant is proportional to the Hubble 
parameter, A{t) oc H{t). We will address these considerations in a future work. 

We also analyse the equivalence between the CV formula and the Friedmann equations with inhomo- 
geneous EoS. Note that this has been done in several works with the homogeneous EoS and the Friedmann 
equations coincide with CV formula only in radiative universe. In this work, we use the inhomogeneous 
EoS including the variable cosmological constant proportional to the square of the Hubble parameter. We 
find that the equivalence between the Friedmann equations occurs only at the present time in a special 
case which is not the radiative universe as for the homogeneous EoS. 

Acknowledgement: The author thanks professors S. D. Odintsov, S. Carneiro and O. Piattella for 
criticism and comments, and CNPq (Brazil) for partial financial support. 
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Figure 1: The Hubble parameter as function of the cosmic time with B = I, from the left to right, A = C = I and 
A = C = —1 respectively. 
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